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Abstract 

We outline a simple way to include heavy quark vacuum polarization 
in lattice QCD simulations. The method, based on effective field theory, 
requires only a trivial modification of the gluon action and has no impact 
on simulation times. We assess the range of validity for this procedure, 
and the impact that it may have. 

1 Introduction 

The inclusion of light quark vacuum polarization in lattice QCD Monte Carlo 
simulations has long been a major problem. With the development of improved 
staggered fermions P it has become possible to do high-precision lattice QCD, 
with good control over all systematic errors ^7. One such systematic error is 
the continued exclusion of heavy quark vacuum polarization. State of the art 
lattice simulations only include the effects of light {s,u and d) quark vacuum 
polarization. 

The exclusion of heavy quark vacuum polarization is a good approximation 
because, as we shall discuss, their effect on simulation results is suppressed by 
a factor of g^/m^, where g is a typical lattice momentum. For many lattice 
calculations (hadron spectra, decay constants) q w ^qcd ~ 250MeV. For a 
charm quark rric « IGeV this could be as large as a 6% effect. Todays lattice 
calculations aspire to few percent (1 — 3%) errors so for truly high-precision 
work, this error should be removed. It would be good to have a method for 
removing it without having to do full dynamical simulations of the heavy quarks. 
Fortunately these effects can be handled using effective field theory. 



2 The Effective Lagrangian 

The earliest example of an effective field theory is the Euler-Heisenberg non- 
linear photon theory 131 . This is an expansion of QED for low-energy photons 
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with average energies uj <C m,, . Euler and Heisenberg were concerned with non- 
linear interactions between four photons, these corrections begin at 0{u}'^/m'^), 
however there is a correction to the photon propagator which is a 0{uj'^ /ml) 
effect 0]. Similar expansions apply to QCD. 

It is straightforward to write down the most general effective action that will 
account for these effects 



S 



EH 



a ^ ^ \ n+1 



On 



(1) 



where C„ are operators of dimension 2{n+l) built out of F^^, F^^ and that 
respect the symmetries of QCD, Cq^ are the coefficients of these operators and 
the sum runs over all heavy quark flavors Q = c,b, t. The operators of higher 
dimensions are suppressed by powers of the heavy quark masses. In this note, 
we will be concerned with only the operators of dimension six. 

At dimension six there are two operators we need to consider |S] 



Gi = gsfabcF^,Fi;,F",^ 
G2 = ^D^FrO^F^, 

With these operators the effective action is 



Seh = I 



(2) 
(3) 



(4) 



This action can be matched to full QCD order by order in perturbation theory, 
since a heavy quark loop is highly virtual. Matching at leading order gives the 
coefffcients [H] 



Q2 



(5) 



3607r' " 157r' 

where Cp = 4/3 for S'[/(3). 

For matching to the lattice theory it is convenient to rewrite the action in 
terms of a different basis of operators. We introduce the operator 



^3 = D^Fl^^'Dy^Fl 



and make use of the relation |S] 
to rewrite the action as 



d'x 



\g^ 



2G2 



Seh — 



d^x 




(6) 



(7) 



(8) 
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In this brief note we are primarily concerned with the effect of heavy quark 
vacuum polarization on the gluon action. With this in mind we make a field 
transformation 

^l^^l + ^^sD-^K, (9) 

which gives 

F^,.Fr - F^^Fr - ^a^g, + O . (10) 

We can pick the constant C to eliminate Q2 from the gluon action^. This 
transformation will introduce additional four-quark operators which should be 
incorporated into the ferniion actions used in the simulations. 

Finally, we rescale the fields A A/gs, and express the action as 

Seh = -ir^ [d*x\Y, Tr [F^.F^,] 

Q Q \ iiu\ / ) 

Our goal is to construct a lattice gluon action that reproduces this action. Before 
wc do this, it is useful to consider how large the mass mg should be in order 
that the expansion in tuq^ is sensible. 



3 Limits on mg 

In order to set limits on the expansion in we look at one part of the 

matching calculation for Q). To match the coefficient Cg^ one can consider 
(light) quark-quark scattering at one loop order in continuum QCD. The only 
effect of virtual heavy quarks at this order will be their contribution to the 
gluon vacuum polarization. Therefore, we can get a check on the validity of the 
expansion by investigating the one-loop gluon propagator in continuum QCD. 
In full QCD the renormalized gluon propagator is given (in Feynman gauge) 

by 

D,Aq') - r'^": V (12) 
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l-n(q2) 

We separate the contribution of heavy quark loops out of the self energy 

n(q2) = ^nQ(g2)^nQCD(g'). (i3) 



^This technique is used to eliminate the same operator in the lattice theory, with m ^ 

„2 
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The light quark and gluonic parts of the self energy will be described by the stan- 
dard QCD action, so we'll just retain the heavy quark contribution. Therefore, 
we will consider the propagator 



(14) 



The contribution to the self energy can be found in any textbook (for example 
0) it's just the QED contribution times the color factor Cp- 



flQ(q^) = -CF— I dxx{l-x)\oi 



1_ J^a;(l - a;) 



Expanding in q^/mg we find 



nQ(<z' 



as q 



(15) 



(16) 



The range of validity for the expansion in 1/toq is set by the radius of 
convergence of the logarithim in (|15|l . 

\q\ < 2mQ. (17) 

Since our goal is to use this expansion in lattice simulations it is useful to express 
this in lattice units, multiplying both sides by a, and taking the momentum 

\q\ = J^QCD 

mqa > O.SAgcna (18) 

For a 0.1 fm lattice h-qcoa ~ 0.1 — 0.3, so we get a very loose limit, mgo > 
0.05 — 0.15. At this lattice spacing, = 0.5, so this bound is satisfied. Notice 
that for low momentum processes it is not necessary for mqa to be greater than 
one, fundamentally this is an expansion in \q\/mQ not l/{mQa). For processes 
that involve larger internal momenta this limit is more strict. For example, the 
semileptonic form factor for B — > tt + ^ + can involve internal momenta as 
high as 



^Qcomb 



IGeV. 



(19) 



With |g| — 1 GeV (|17|l is still satisfied for a charm quark, but for lighter quarks 
it wouldn't be. When we match to the lattice theory we will find another bound 
on mna. 



4 Matching to the Lattice Theory 

The leading correction in (|ll|l is trivially included in the widely used improved 
gluon action 

Sl = ^Y1 |/3p^ReTr [1 - Up] + /SR^RcTr [1 - Ur] + /?6^ReTr [1 - Ug] 



(20) 
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where Up is the plaquette, Ur is the rectangle, and Uq is the parallelogram 
term^. We define the operators 

To = ^Tr(F,,,F,,,) 

= ^ Tr (D^F^^D^F,,,) (21) 

in terms of which we have 

iReTr[l-Up] = -^T, + ^^T, + O {a^) 

iReTr[l-UR] = _2a4^o + + O (22) 

iReTr[l-U6] = -2a^To +'^^2 + O {a^) . 



The choice^ 



/36 = 

Pn = (23) 



1 



gives [TU] 

5^ = 



--^ / d^xTo + O (a^a', a^ ^^/(m.a)^) . (24) 
2.90 J 

In order to reproduce (fTT|l we need to have a non-zero value for /Jg- When 
we do this, we are be forced to included an additional correction to [3r in order 
to cancel any potential contribution from Ti. We take 

ul [arnqY 



■^Corrections to the standard Wilson plaquette action due to heavy quark vacuum polar- 
ization have been investigated in Isllol . 

''The factor of uq is a tadpole improvement factor, it plays no role in our discussion. 
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which gives 



250 V 



1 - 



{amqf 



15cw - ^Rl'l 

3 3 



^1 



9(aTOQ)2 

where for simplicity we just include only one flavor of heavy quark. Setting 

1 1 



(26) 



^[1] 



2007r 



fSOTT 



(27) 



produces the correct coefficient for the ^2 term. The additional contribution to 
the coefficient of the J-q term produces an unobservable shift in the wavefunction 
renormalization for the gluons, so we are free to drop it. This leaves only the 
J-i term. We tune to insure that its coefficient remains zero. This is easily 
accomplished with 

= f2C7[il = - (28) 

ZOTT 

With these values the lattice action correctly reproduces Ullll with 
O (a"*, a^a^, tts/ (arric)'^, al/{amcY) 



(29) 



errors. 



5 Impact on One-Loop Corrections 

We combine our results with the known one loop corrections to Pr and (3q, 



I3f 



20ul 



I + as (TT/a) I 0.4805 + XnNf - ^ E 



257r ^ (amgy 



(30) 



/3e = -|..(V«)(^0.03325 + XeiV, + ^E(^j- (31) 

The one loop terms due to the gluons are well known ^2 ^| j those due to dy- 
namical light fermions are being computed jl3|. With amc ~ 1 these corrections 
to the one-loop terms are 0.025 for P^ and 0.0021 for P^. 

These expressions also limit how small one can take a. The coupling constant 
as {t^ / a) goes to zero logarithmically as a — > 0, whereas the correction from the 
heavy quark vacuum polarization is growing as Clearly, if we take a 

too small the correction term will end up larger than the tree level term. For 
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example, on an ultra-fine lattice, with a = 0.01 fm, we have rtica = 0.05, which 
gives (setting Nf = 0) 



/3fl = -:^[l-9-7a«(7r/a)], /36 = -^G.88a, (7r/a) . (32) 

In this case, the "correction" to I3r is larger than the leading order term. With 
TOcfl = 0.13 the coefficient of the correction term becomes one, so we take this 
as a rough lower limit, a ~ 0.025 fm. 

A more realistic lattice spacing is a = 0.09 fm, which gives rrica = 0.45 and 



I^R = + (0.481 -0.126)a,(7r/a)] = --^[l + 0.355a47r/a)], 



20u^ ' ' / - V / /J 20u^ 

I3q = -^(0.0333 + 0.0105)as (7r/a) = -^0.0437as(7r/a). (33) 

In the this case, the corrections due to heavy quark vacuum polarization are 
significant but do not cause convergence problems. 



6 Conclusion 

In this note we have demonstrated that heavy quark vacuum polarization can be 
included in Monte Carlo lattice simulations with no additional computational 
cost. The action is the standard Symanzik improved gluon action 

Sl = ^Y1 [(^pI^^^' [1 - Up] + /SR^ReTr [1 - Ur] + /Se^ReTr [1 - U6]| 

(34) 

with modified coefficients 

in /a) (|o.4805 + X,Ar,-J-X:^j 
0.03325 + XeiV, + ^X:(^) (3^) 

This action reproduces continuum QCD up to corrections of order a'',(asa)^, 
as/(amc)^ and al/{amc)^. In order to insure reasonable behaviour of the per- 
turbation series rUca > 0.13 is required. With this constraint, this action should 
simulate processes with average internal momentum |g| <C 2mc. Ideally, the 
charm quark would be treated dynamically, then these constraints would apply 
ruba instead. 



Pr = 
= 

pp = 



20m§ 



1 + Qfs 



— jas (TT/a) 
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